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ABSTRACT 

We present a detailed analysis of a disc galaxy forming in a high-resolution fully cosmo- 
logical simulation to investigate the nature of the outer regions of discs and their relevance 
for the disc formation process. Specifically, we focus on the phenomenon of misaligned disc 
components and find that the outer disc warp is a consequence of the misalignment between 
the inner disc and the surrounding hot gaseous halo. As the infalling cold gas sinks toward 
the centre of the galaxy, it is strongly torqued by the hot gas halo. By the time the fresh gas 
reaches the central disc-forming region its angular momentum is completely aligned with the 
spin of the hot gas halo. If the spin of the hot gas halo, in turn, is not aligned with that of 
the inner disc, a misaligned outer disc forms comprised of newly accreted material. The in- 
ner and outer components are misaligned with each other because they respond differently to 
infalling substructure and accretion. The warped disc feeds the main gas disc due to viscous 
angular momentum losses, but small amounts of star formation in the warp itself form a low- 
metallicity thick disc. We show that observations of resolved stellar populations in warped 
galaxies in the local universe could provide evidence for the presence of these processes and 
therefore indirectly reveal ongoing gas accretion and the existence of hot gas halos. 

Key words: galaxies: evolution — galaxies: spiral — galaxies: stellar content — Galaxy: 
solar neighborhood — Galaxy: stellar content — stellar dynamics 



1 INTRODUCTION 

A large fraction of the baryonic mass in a spiral galaxy resides in 
the disc component. The discs of massive spirals assemble pri- 
marily from the quies cent cooling of shock-heated gas after the 
last major merger (e.g . iFall & Efstathioijl 19801 : ISrook et al.ll2004 
[Robertson et al. l2006h . although depending on the mass of the 
galaxy, cold accretion also plays an important role in providing 
disc gas (e.g. Keres et al. 2005; Brooks et al. 2009). Lower angu- 
lar momentum material is initially denser and closer to the cen- 
tre of the galaxy, therefore cooling faster and accreting early. It 
may also be generate d through angu l ar mo mentum loss during 
the merging process jP'Onghia et"^l2006h . which in the hier- 
archical structure growth paradigm is more vigorous during the 
early stages of mass assembly, and yields old spheroids in the 
centers of galaxies. Higher angular momentum material follows 
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late, resulting in inside-out growth of galactic discs ( lLarsod[T97^: 
|white & Frenk 1991). Observationally, inside-out growth can be 
inferred fr om trends of d e creasing mean stellar age as a function 
of radius Me long" 1996'; 'Bell & de Jon3 1200(]|; iMacArthur et al] 



^004; Williams et al. 2009; Go garten et al 
ment of specific star formation rates as 



2010l) . or by measure- 
a function of radius 



Muiioz-Mateos et al. I l2007l) . By extension, the outermost regions 
of discs are the sites of current disc assembly, and therefore invalu- 
able laboratories for studies of spiral galaxy formation. 

The outer reaches of discs are inherently difficult to study ob- 
servationally due to low surface brightnesses several disc scale- 
lengths from the galactic centre. On the other hand, the high numer- 
ical resolution needed to study their properties hampers theoretical 
studies. However, it has long been acknowledged that these extreme 
regi ons of galaxies m ay hold unique clues to galactic disc forma- 
tion. |Freem^ {1973) realized that the surface brightness profiles of 
some discs defy the usual ( and simple!) characterization by a sin- 
gle exponential jde Vaucou leurs 1958) in their outer parts. Subse- 
quent observations have revealed that downward-bending exponen- 
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tial profiles are very common, while pure exponential and upward- 
bending profiles constitute a srnaller fraction Ivan der K ruilll 19791 : 
iKreeel et alj|2002l : IPolilen et alJ20o3: IPohlen & Truiillol20 0#). 

The elusive outer discs present another puzzle - the extended 
HI is rarely aligned to the same plane as the inner disc. Misalign- 
ments between different components of spiral galaxies have been 
studied exte nsively over the past several decades (Briggs 1993; 
lRubinll994h . Phenomenologicall y, such misalignmen ts range from 
warps, usually observed in HI i Sancisil 1 19761 : iGarcia-Ruiz et al., 
l2002l:IVerdes-Montenegro et alj|2002h " to structures whose angular 
momentum axes are orthogonal to that of the main disc, i.e. "polar 
disc/ring" galaxies (e.g. Whitm ore et al , 1990). In a modest sample 
of inclined galaxies. iGarcfa-Ruiz et aLr0OO2h found that all galax- 
ies whose HI discs extend beyond the main stellar disc also have a 
warp. Warps have been shown to origina te through torquing of the 
disc either by misaligned cosmic infall tostriker & Binney" 19891; 
guinn&Binnev 1992; Jiang & Binnev 1999; Shen & Sellwoi3 
20061) . misalignment betw een the angular momentum of the disc 



and the dark mat ter halo jPebattista & SellwoodI 19991) . or satel- 
lite perturbations dWeinbe rg & Blitz 2006). The origin of polar ring 
structures has been the subject of much debate and there are many 
plausible mechanisms for their formation. Polar discs are often ob- 
served around SO sp irals and are postulated to form d uring an inter- 
action I Bekkjl 1998h or cosmological gas accretion jMaccio et al.l 
l2006l ; lBrooketalj2008l) . 

Comparatively little work h as been done on the stellar c onten t 
coinciding with misaligned HI. IVerdes-Montenegro et al. I j2002h 
showed that the warped outer layers also support low-level star 
formation in the case of NGC3642. Recent results from GALEX 
show faint UV emissi on in the warped discs of M83 ( Thilker et al. 
12005 ) and NGC5055 jThilker et al.ll2007l ; ISancis"ieralJl2008l) . im- 
plying that a low level of star-formation also exists in those warps. 
In th ese cases, the UV emission is found well beyond the Ha "cut- 
off" (lKennicuttlll98^ . resulting in an extended-UV (XUV) disc 
and at least circumstantially linking HI warps to this phenomenon. 
Systematic studies of optically-identified warps agree with HI es- 
timates of warp frequency and suggest that > 50% of all galaxies 
are w arped jSanchez-Saavedra et alJl990l ; lReshetnikov & Combed 
Il998h . 

The warps discussed above can only be observed in the local 
universe, but their ubiquity suggests that they are long-lived phe- 
nomena (unless they are all caused by recent satellite interactions). 
If large-scale HI warps harbor sufficiently dense gas to allow for 
star formation, several generations of stars must have been born out 
of the plane of the main star forming disc. Furthermore, since discs 
grow with time, stars bom out of the plane of the main disc must 
contribute to the stellar population not only in the tenuous outer re- 
gions, but throughout the disc, possibly contributing to a thick disc 
population. We address these intriguing possibilities below. 

We have previously addressed oute r discs in the context of disc 
dynamics and stellar radial migration in jRoskar et alj2008al lbl R08 
collectively hereafter). In this Paper, we focus on the importance of 
cosmological gas accretion and the intrinsic misalignments of dif- 
ferent components for the outer parts of galactic discs. We study in 
detail the underlying cause for the misalignment between the cen- 
tral disc and the newly accreted gas, and find that it is the overall 
angular momentum of the entire gas reservoir which enables the 
misalignment. We therefore propose another explanation for the 
ubiquity of warps, and link this common feature of disc galaxies 
to their broader cosmological context. 

The Paper is organized as follows: in §|2]we briefly describe 
the simulation and the relevant details of the code; in §[3] we de- 



scribe the morphological properties of the misaligned disc; in §|4] 
we describe the physical origin of the misaligned disc in the sim- 
ulation; in §[5] we describe the properties of the stars that form in 
the misaligned disc; in § |6] we discuss the possible influence of 
long-lived warps on thick discs; in §|7]we show that the processes 
described in the preceding sections may be detectable by resolved 
star studies of nearby systems; in §[8]we discuss potential overlap of 
the warp signature imparted on the stellar populations to that of ra- 
dial migration, and show that the disc is too hot to harbor sufficient 
non-axisymmetric structure needed for such secular evolution; we 
summarize our findings in §|9] 



2 METHODS 

The simulations presented in this Paper were run using the A'^-body 
+ SPH parallel tree code GASOLINE (Wadsley et al. 2004), using 
Compton and radiat ive cooling, star form ation and supernova feed- 
back as described in'Stinson et alj ( l2006h . and a UV radiative back- 
ground jHaardt & Madau 1996). The star formation recipe follows 
the prescription from lKatd ( I1993) . forming stars according to 



dp. 



idyn 



where we use c* — 0.05. pgas and tdyn are the local gas density 
and dynamical time evaluated using the SPH kernel. A gas particle 
is eligible to form stars if its density exceeds 0.1 amu/cm"^ and its 
temperature is below 1.5 x 10* K. The feedback recipe accounts 
for both SN II and SN la explosions, keeping track of iron and oxy- 
gen abundances. Supernova explosions impa rt energy u pori p arti- 
cles enclosed within a blast radius based on iChevalieil ( ll974l) . To 
prevent the particles from immediately radiating away the energy, 
radiative cooling is turned off for a time corresponding to the end 
of the snow-plow phase of the SN remnant. 

We use two simulations in this paper. Our fiducial simu- 
lation, used to describe most of the physical proc esses in jj [3l 



Brooks et al 



Brooks et al 



§ [4] an d [5] is the MWl run discussed a lr eady in 
120071) : iPontzen et all j2008l) : iRead et al.1 ( l2009h ; 

(2009). A lower-resolution version of the simulation was presented 
in Gover nato et al. (2007, 2008). The reader is referred to those 
works for details on the initial conditions, here we just briefly out- 
line the salient points. The simulation was run with WMAP year 
1 cosmological parameters for the A-dominated, flat cosmology: 
Qo = 0.3, A = 0.7, ilt = 0.039, h = 0.7. At z=0, the total mass 
within R^ir — 270 kpc is Mvir = 1.1 x 10^^ Mq, resolved by a 
total of 1.26 X 10'^ , 2.85 x lO" and 6.7 x 10^ dark matter, star, and 
gas particles respectively. The particle masses are 1.4 X 10^ Mq 
and 7.6 x 10^ Mq for the gas and dark matter respectively. Stars 
form with a mass of 4.8 X 10'' Mq. The softening lengths of all 
particles are 0.3 kpc. Low-temperature cooling processes were not 
included in this simulation, limiting the minimum gas temperature 
to ~ 10* K. 

Our fiducial simulation described above is an example of a 
system that was warped in the past but is entirely co-planar at z = 
0. Therefore, we focus in §|7]on another simulation (h603), which 
preserves its large-scale warp to the present day. The simulation 
was also run with GASOLINE, using WMAP year 3 cosmologi- 
cal parameters (Qo = 0.24, A = 0.76, h = 0.73). Its total mass 
within Ryir = 186 kpc is Myir = 3.7 x lO" Mq. The initial con- 
ditions were sampled on a 2304^ grid, resulting in 5.4 x 10^ gas, 
2.2 X 10® star, and 1.0 x 10® dark matter particles within the cen- 
tral galaxy's virial radius at z = 0. Initial particle masses for the gas 
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and dark matter particles (in the central region) are 6.3 X W Mo 
and 3.0 x 10^ respectively. Initial stellar masses are 2.1 x 10^ 
Mq. Softening lengths are 0.23 kpc for all particles in the inner re- 
gion. This sim ulation uses a comb ination of low-temperature metal 
cooling ( Mash chenko et al1l2006h and self-shielding sub-grid pre- 
scriptions, both of which are mostly relevant for evolution at high 
redshift, which is not the focus of this Paper. The star formation 
recipe and its parameters are identical to the fiducial MWl run. All 
sections that follow with the exceptions of §|7]and §[9]refer only to 
the fiducial run, MWl. 



t = 5.76 Gyr 



t = 6.90 Gyr 





t = 10.10 Gyr 



t = 13.47 Gyr 



3 MORPHOLOGY OF THE MISALIGNED DISC 

A misaligned disc is evident in our fiducial simulation almost im- 
mediately upon the completion of the last major merger at z ~ 2, 
and persists for ~ 7 Gyr. In Figure [T] we show HI column density 
maps showing the central disc and the extended warp at z 0.5. 
To measure the tilt, we rotate the simulation so that the angular 
momentum axis of the gas in the inner 3 comoving kpc is the z- 
axis. We then measure the angular momentum of gas in concentric 
annuli and take the tilt to be the angle between the angular momen- 
tum axis of the annulus and the jz-axis. Note that this definition of 
warping differs from some de finitions used in the litera ture to char- 
acterize observed warps (e.g. iGarcia-Ruiz et alj|2002h . but we are 
not seeking exact model matches to the data; the measurement of 
tilt is purely for the determination of misalignment in the simula- 
tion and defining the tilt in terms of the angular momentum axis is 
a natural choice. In the analysis that follows, we define the warp to 
begin at a radius where the tilt angle exceeds 10°. 

In Figure |2] we show the orientation of the angular momen- 
tum vectors of concentric radial bins with respect to the main star- 
fo rming d i sk us ing Briggs figures. Briggs figures were introduced 
bv lBriggsl i fT99Q) to easily visualize the radial evolution of a warp 
by representing each radial bin as a point whose radial coordinate 
is the tilt (relative to the angular momentum of the central region) 
and the polar coordinate is the angle of the line-of-nodes (LON) 
from the a;-axis. The LON is the imaginary line where the warp in- 
tersects the ly-plan^B- From this representation, it is clear that the 
misalignment can in fact be considered a warp, given that succes- 
sive rings twist and tilt with increasing radius. If, for example, the 
outer misaligned disk was an independent flat structure, we would 
see two clusters of concentric points in Fig. [J] The warp is already 
quite pronounced at ~ 5.7 Gyr, reaching a tilt angle of ~ 40° by 
7 Gyr. The warp gradually subsides and almost completely disap- 
pears by the present-day (we discuss the reasons for the disappear- 
ance of the warp in §|43J. When a strong warp is present, it always 
traces out a leading spiral in the tilt-L ON representa tion, consistent 
with the properties of observed warps terigg^ 199(f). 

Our warp definition necessarily fails to identify a weak warp 
(such as in the bottom right panel of Figure |2}, but we find it im- 
possible to reliably separate the warped layers from the un warped 
disc when the tilting is < 5°. A perfect decomposition is impossi- 
ble because the transition between the two components is smooth, 
and the misaligned disc in fact feeds the main star-forming disc. 
However, as we show below, the gross properties of the large-scale 
warp that persists for a substantial portion of the simulation can be 



^ The polar angle shown is actually the angle between the angular momen- 
tum vector of the bin and the x-axis, which is 90° from the LON. However, 
since we are primarily interested in the change from one radius to the next 
rather than absolute angles, this difference in definition is irrelevant. 



Figure 2. Briggs figures at several different times. Each point represents a 
radial bin. There are a total of 15 bins spaced equally between and 10 kpc. 
The radial coordinate in the polar representation is the angle of the angulai' 
momentum axis of the bin from the angular momentum axis of the main 
star-forming disc. The angular coordinate is the line of nodes (LON). The 
dotted circles ai'e plotted in intervals of 10°. Direction of disc rotation is 
counter-clockwise. 



related to large misalignments between different galactic compo- 
nents, the understanding of which is not predicated upon a precise 
distinction of warped and un warped parts of the disc. 



4 ORIGIN OF THE MISALIGNED DISC 

Understanding the physical cause of the misalignment is crucial be- 
cause of its potential to link an observational property of a galaxy to 
its broader, cosmologically relevant surroundings. Dynamical stud- 
ies have successfully produced warps in simulations, by subject- 
ing a flat disc of concentric rings or particl es to an ext ernal torque 
(Ostriker & Binnev 1989; Ouinn & Binne^J|l99^^^img & Binne; 
1999; Debattista & Sellwood 1999; Lopez-Corredoir a et alWim 
,Shen & Se llwood 2006; gan chez-Salced_Q 2QQ6 ). In the case of our 
simulation, the misaligimient precedes the evolution in the main 
disc plane and is a consequence of the gas accretion. 

Several recent theoretical studies have addressed the process 

c J r ^ 1 ■ F^- — Z~- P 1 I'^nm' li^ ' 



jue 
ie3 
qI; 



2005 


; Dekel & Bimboiml 


2006.; .Keres et al., ,200a; ,Ocvirk et al. 


200S 


;lDekeletal, 2009allb 


; Brooks et al. 20091. hereafter B09). If 



the halo is massive enough (> IO^^Mq) to support a stable shock, 
some gas is shock-heated to the virial temperature upon entering 
the halo (hot accretion), while other gas penetrates the halo in- 
side a filament and is never shock-heated (cold accretion). Because 
such cold flows could potentially deliver substantial amounts of 
specific angular momentum deep into the galactic potential well, 
misaligned components of disc galaxies have been invoke d pre- 
viously as manifestations of cold accretion fcaccio etal] |2006l ; 
iBrook et al.ll2008l ; ISancisi et alj2008l ; ISpavone et alj2010il . 

While most of the above studies use a temperature criterion to 
separate the two modes of smooth accretion, B09 additionally used 
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Figure 1. HI maps of the disc at 2 ~ 0.5 in a few different orientations. Colour indicates logarithmic HI column density, distance scale is in comoving kpc. 
The main disc is seen as the high-density central region in the inner ~ 10 kpc. 
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Figure 3. Total gas mass in the misaligned disc and main star-forming disc, 
separated by mode of accretion. Red (shocked) coiresponds to particles that 
enter the virial halo and shock at the virial radius, blue (unshocked) shows 
particles that enter the virial halo unshocked, and green (clumpy) shows 
particles that enter the virial radius after having belonged to another halo. 



an entropy criterion on the basis that a strong shock will always 
induce an entropy jump, but not every high-temperature particle 
must necessarily have been shocked. In what follows we use the 
same particle groupings as B09, dividing the gas particles in the 
main galaxy halo into unshocked, shocked and clumpy (accreted 
after being a part of another halo). These definitions are applied to 
particles more than 30 kpc from the central star-forming disc, so 
as to not confuse entropy jumps and temperature changes due to 
large-scale accretion shocks with those from supernova feedback 
within the star-forming disc. Below we explore whether the mis- 
aligned gas component in our simulation arises due to cold flows, 
dark matter halo torques, or whether it results from another mech- 
anism. 

In Figure [3] we show the gas mass found in each compo- 
nent (misaligned disc and main star-forming disc) as a function of 
time, broken down according to the mode of accretion, as discussed 



above. The total amount of gas found in the misaligned component 
is as much as a factor ~ 2 larger than that in the main star-forming 
disc. Importantly, though the misaligned disc is mostly composed 
of gas acquired through cold accretion {r^ 75 per cent), all modes of 
accretion contribute to its total mass. The misaligned gas therefore 
cannot only be a consequence of accretion of unshocked gas, but 
must also be due to a more general mechanism that affects all gas 
equally. In the following sections, we explore possible mechanisms 
for this behavior. 



4.1 Warps as Indicators of Misaligned Hot Gaseous Hales 

In this section, we propose that the outer disc warps are a result 
of the misalignment between the inner disc and the surrounding 
hot gas halo. We define "hot" and "cold" gas as all gas with T > 
10^ K and T < 10^ K respectively. As the cold gas sinks toward 
the centre of the galaxy it is torqued by the massive hot gas halo 
long before it reaches the disc. Due to this torquing in the outer 
halo, the cold gas aligns with the spin of the hot gas halo by the time 
it reaches the central region. Because the hot gas halo is misaligned 
with the main star-forming disc, its misalignment manifests as a 
warp that is mainly composed of cold infalling gas. 

In Figure|4]we show explicitly that the hot gas halo is responsi- 
ble for torquing the infalling cold gas. In the top panel, we show the 
specific torque resulting from hydrodynamic forces on cold and hot 
gas (solid and dashed lines respectively). The specific torque on the 
cold gas is much higher in the outer parts of the halo, because the 
hot gas dominates the cold gas by mass in the outer regions. How- 
ever, the total torque on each subset is almost exactly the same, 
as shown in the bottom panel. We have verified that the directions 
of the SPH torques on the two types of gas are antiparallel in the 
outer parts of the halo. Hence, the hot and cold gas exert a torque 
on each other, but the more massive hot gas halo dictates the sense 
of rotation for the infalling cold gas. Although we only show one 
timestep in Fig. U this behavior is evident at other timesteps we 
analyzed after z = 2. 

We have also investigated the evolution of hot gas halo gas 
present within the virial radius as a function of time. Although the 
angular momentum of the hot gas halo gas changes through accre- 
tion, the halo particles behave coherently. For gas that makes up 
the hot gas halo (T > 10^ and 0.3 < r/r^rr < 0.7) at 5 Gyr, the 
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Figure 4. Top and bottom panels show the specific and total hydrodynamic 
torques respectively on cold (solid lines) and hot gas (dashed lines) at z 
~ 0.5 (t ~ 8 Gyr). The specific torques are much different because hot gas 
dominates by mass in the halo, however the absolute torque magnitudes are 
equal, showing that the hot gas torques the cold gas. 



angular momentum direction remains largely radially coherent at 
all times. 

The torquing of cold gas results in very close alignment of an- 
gular momenta of the hot gas halo and infalling cold gas. In the top 
panel of Figure[5]we show the radial profile of 6 (solid lines) and 
(j) (dashed lines), angles between the z and x axes of the main star- 
forming disc and the angular momentum vector of gas in each bin 
in a single snapshot. The gas is subdivided into cold and hot, again 
based on a temperature criterion regardless of whether it shocked 
or not, as we are interested in seeing whether cold gas preserves its 
trajectory through the halc0. The bottom panel explicitly shows the 
alignment of the hot and cold gas defined by 



Jl-J2 

\Jl\\J2\ 



(1) 



The top panel shows much more variation in the profile of the cold 
gas, though on average the cold gas follows the hot gas, especially 
closer to the centre of the halo. The angle between the two compo- 
nents steadily decreases inwards (bottom panel) and they converge 
completely in the inner region where the disc fueling takes place. 

In Figure [6l we show the evolution of particles that enter the 
virial radius at t ~ 3.8 Gyr, grouped according to their angular 

^ Selecting particles purely based on mode of accretion may confuse the 
analysis because it could include particles accreted as unshocked, but ex- 
pelled from the disc through supernova feedback and therefore belonging 
to the hot gas halo. By limiting the selection according to present tempera- 
ture, we make sure in part that the cold particles are ones currently falling 
toward the central disc for the first time. These may and do include particles 
that shocked upon accretion. 
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Figure 5. Top: Mean angles 8 (solid lines) and </) (dashed lines) of gas an- 
gular momentum vector subdivided by gas temperature. Note that although 
the gas with lower temperature will sink quicker to the centre of the halo, 
the hot gas dominates by mass and drives the sense of rotation. The region 
of the misaligned (warped) disc is roughly between 0.05 < r/ryir < 0.15. 
Bottom: Angle between the cold and hot gas angular momenta as a function 
of radius. 



momentum aligimient upon entering the halo. We chose two very 
different groups of particles in order to investigate whether the in- 
coming angular momentum is preserved as the gas sinks to the cen- 
tre. We follow the mean properties of the particles in each timestep, 
though the interpretation is made more difficult by the fact that the 
particles do not necessarily remain grouped spatially. Nevertheless, 
this approach sheds some light on the gas evolution within the virial 
radius. 

Specifically, we select the particles according to their angular 
momentum vector with respect to the angular momentum of the 
central disc, yielding two subgroups related both spatially and by 
their direction of motion. Top and bottom panels show particles 
that were co- and counter-rotating as they entered the virial radius 
respectively. The range of angular momentum vector angles used 
for the selection are 118° < < 122° and -160° < (j)< -140° 
for counter-rotating particles and 70° < 8 < 78° and 50° < (f) < 
70° for the particles co-rotating with the disc. 

Even in the extreme case of unshocked counter-rotating gas 
entering the halo (bottom panels, blue lines) the sense of rotation 
with respect to the central disc is not preserved, indicating that the 
infalling gas experiences substantial torques as it travels through 
the halo. To demonstrate that the hot gas halo is indeed the source 
of the torques, we compute the rotation of the hot gaseous halo 
within a shell of thickness 0.05 r^ir about the mean radius of the se- 
lected particles, and show it in the right-hand panels with diamond 
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Figure 6. Time evolution of particle properties for gas particles that enter 
the virial radius at the same time (t ~ 3.8 Gyr). Top and bottom pan- 
els show particles that upon entering the halo are co- and counter-rotating 
respectively. The particles were chosen from narrow regions of 9 and (j) dis- 
tributions, where and specify the angular momentum vector of the in- 
dividual particles with respect to the central disc (see text for more details). 
Black, blue, and red lines represent all, unshocked, and shocked particles 
respectively. Symbols in the right panels denote the rotation of the halo 
around the mean radius of the particles at a given time. Both sets of par- 
ticles, co-rotating and counter-rotating, eventually align with the angular 
momentum of the hot gas halo. 



symbols0 The infalling cold particles in general follow the halo 
rotation quite closely, in agreement with the expectations based on 
Fig.g] 

Figure [6] also shows that the gas that entered the halo in the 
hot phase sinks to the centre more slowly, and a smaller fraction 
makes it to the star forming region (left-hand panels of Fig. 
About 50 per cent of the unshocked gas and about 30 per cent of the 
shocked gas gets consumed by star formation before being heated 
by supernova feedback and ejected back into the hot gas halo. Once 
the gas is ejected from the disc, it becomes a part of the hot gas halo 
at temperatures > 10^ K, which is evident by the rapid increase in 
mean r. Also note that the large discrepancy in the alignment at 
8-9 Gyr is due to an infalling satellite which complicates a robust 
determination of angular momentum angles at those times. 

We note that a qualitati vely similar phenomenon was pointed 
out by iKatz & Gunnl ( Il99 ih . who also found that the outer parts 
of the gas disc were warped with respect to the centre, but aligned 
with the angular momentum vector of the initial top-hat distribu- 
tion. They attributed the offset to the differences in the details of an- 
gular momentum transfer between the gas blobs and the dark mat- 
ter halo. The "angular momentum catastrophe" in today's state-of- 
the-art simulations has been alleviated by increased resolution and 
more sophisticated modeling, but in the fuUy-cosmological context. 



This is not ideal because, as we noted above, the particles we select ini- 
tially do not necessarily remain tightly grouped. For example, the offset in 
6 between the shocked and unshocked particles at ~ 8 Gyr is due to them 
occupying slightly different parts of the halo. Nevertheless, this comparison 
is useful for relating the infalling particles to the halo in a general sense. 



the added complications of dark matter halo triaxiality and merging 
facilitate the perpetuation of misalignment. 



4.2 Misalignment due to Dark Matter Halo Torques 

Though we have shown above that the hot gaseous halo torques 
the infalling cold gas, the torquing by the gravitationally dominant 
DM halo could dictate the overall angular momentum properties of 
the gas and ultimately, the misaligned disc. Two main possibilities 
exist: a) the DM efficiently torques the gas particles, influencing 
their angular momentum alignment in the halo and b) the DM halo 
torques the flat disc into a warp. We address each of these in turn 
below. 

Investigating the first of these two possibilities, we first mea- 
sure the shape of the dark matter halo in ellipsoidal shells using an 
iterative procedure. On the first iteration, the mass moment tensor, 
defined by 



(2) 



in a given spherical shell is diagonalized to obtain the principal axis 
vectors and the associated eigenvalues. Note that this is not the mo- 
ment of inertia te nsor, although in the literature this nomenclature 
is often used (see iBett et all 2007h . On subsequent iterations, we 
redefine the ellipsoidal shell according to the eigenvectors by 



rell 



y 

b/a 



+ 



c/a 



(3) 



after rotating the particle coordinates into the frame where the mass 
distribution matrix is diagonal. The axis lengths a, b, and c are 
given by the square roots of the eigenvalues. Once the shell is re- 
defined, the mass distribution matrix is reconstructed and the pro- 
cedure repeated. We find that in most cases the deviations in axis 
ratios are less than 1 per cen t after ^ 2 iterations. This procedure 
is analogous to that used in rKal3 h99lh. though it uses differential 
instead of cumulative shells following Debattista et alJ l l2008 h. 

For the majority of the simulation from 2 r^ 2 to the present, 
the axis ratios of the DM halo are b/a ^ 0.8 — 0.9 and c/a ~ 
0.5— 0.7 beyond r /r^ir ~ 0.2. We find that in the innermost region 
the angular momentum of the gas is well-aligned wi th the minor 
axis of the DM halo (also reported by Bailin et al. 2005), but we do 
not find any obvious coupling between the gas angular momentum 
vector and the dark matter halo's minor axis during periods when 
the warp is present. 

To check explicitly the importance of DM torques on the tra- 
jectories of gas particles through the halo, we calculate gravita- 
tional and hydrodynamic (SPH) torques on the particles directly. 
Figure |7] shows the comparison of torques from DM and hydrody- 
namics on the cold and hot gas. Recall that the total SPH torques 
were equal on these two subsets of particles (Fig.|4}. Fig.|7]shows 
explicitly that the DM torques are insignificant in shaping the tra- 
jectory of gas particles through the halo. Differences between the 
hot and cold gas arise because the cold gas is in general patchy 
and irregularly spaced in each radial shell, therefore experiencing 
a larger net gravitational torque in each shell from the dark matter 
halo than the more homogeneous and spherical hot halo gas. 

This result seems counter intuitive, given that the dark matter 
dominates the potential of the galaxy, but a simple experiment con- 
firms that such behavior is expected. To simplify the problem, we 
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Figure 7. The ratio of SPH torque to DM torque on the cold and hot gas in 
logarithmically-spaced radial bins at z~0.5. Except for the innermost bins 
in the region of the central disc, the SPH torque dominates over that of the 
dark matter. 



imagine a triaxial halo of uniform density at the origin, with princi- 
pal axis lengths a, b, and c (aligned along x, y, and 2 axes respec- 
tively), where a > b > c. We assign a mass Mhaio = 3 x 10^^ Mq 
to the halo and set the longest axis to be a = Q.^Vvir- The other 
two axis lengths are set according to axis ratios measured in the 
simulation with b/a = 0.8 and c/a = 0.6. We then place rings 
of test masses at a distance Tc to represent the clouds of infalling 
material. The two rings are inclined with respect to the xy-plane 
by 45°, rotated around the x and y axes respectively. We may then 
calculate the gravitational torques Tg on the test masses in the rings 
using equations for potentials of homogeneous spheroids given in 
Table 2. 1 and 2.2 of Binney & Tremaine ( 2008.) . 

Similarly, we can make a simple estimate for the expected ram 
pressure force on a cloud of gas moving through the rotating hot gas 
halo. The torque due to the ram pressure Pram ~ phV^ei is given 



by 



(4) 



We assume in this estimate that the clouds in our simulation are 
analogous to compact high velocity clouds (CHVCs) observed 
in the halo of the Milky Way and other nearby systems (e.g. 
IWestmeier et al Jl 20051 : IPutman et alJl2002l : |Pisano et al.ll2007l) . We 
assume that the clouds are moving with velocities 100-200 km s^^ 
relative to the hot gas halo. Sizes of CHVCs are fairly uncon- 
strained due to inherent difficulties in estimating their distances, 
though most constra ints place HV Cs within 150 kpc of the host 
galaxy centers (see IPutmarj 1200^ . To estimate the si ze of the 
clouds , we take the mean angular size CHVCs from Put man et al.l 
Asky = 0-36 deg'^. Assuming that the clouds are at 150 
kpc, we obtain a radius of ~ 1 kpc. We consider this to be an up- 
per limit and consider radii in the range of 0.05 - 1.0 kpc. Using 
equation |4] and placing the clouds anywhere between O.Gvvir < 
< Tvir, we find that Tg/rram ^ 30 per cent for reasonable 
choices of halo and cloud parameters. The behavior seen in the 
simulation, where gas trajectories are determined largely by hydro- 
dynamic torques, is therefore expected for cold gas encountering a 
hot gas halo in the outer parts of the galaxy. 

We now try to determine whether the warp could be caused 
by the torques of the outer halo on a flat disc, as was re- 
cently demonstrated as a plausi ble warp formation mechanism by 
IPubinski & Chakrabartvl ( l2009l) . We calculated the torque in the 
disc plane due to the DM halo particles beyond r/rvir > 0.2 (i.e. 



beyond where the halo and the disc are aligned) in a few timesteps 
where the warp is clearly evident, including the output shown in 
Figure [T] The timescale for changing the angular momentum J of 
a ring at radius r by A J ~ J is given by t{r) — J{r) /t±, where 
J is the angular momentum of the ring and t± is the component 
of the torque perpendicular to the ring's angular momentum vector. 
We find that this timescale is > 10 Gyr everywhere. Considering 
that we see a tilt of ~ 40°, J only needs to change by 25 per cent, 
but the timescale for such change is still several Gyr in the region 
of the warp. If this were the dominant mechanism for the warping, 
we should first see a flat disc develop, whose outer regions would 
torque out of the plane slowly. This is not the case, as the warp is 
in place almost immediately after the last major merger. Therefore, 
the outer DM halo torque cannot produce the warp in this system. 

We conclude that warp formation in this simulation is not sig- 
nificantly affected by the dark matter torques because hydrody- 
namic torques dominate in the outer halo where the spin of the 
infalling gas is determined prior to its entry into the disc region. By 
the time the fresh gas reaches the central region, it is aligned with 
the hot gas halo, but misaligned with the central disc, giving rise 
to a warp. The warp is therefore a manifestation of a misalignment 
between the spin axes of the central disc and the surrounding hot 
gaseous halo. 



4.3 Onset and End of the Warp 

We have established in the previous section that the misaligned cold 
gas component in our simulation is due to distinct spin axes of the 
inner disc and the surrounding hot gas halo. Here, we try to deter- 
mine what causes the two components to be misaligned in the first 
place and why the warp eventually disappears. 

First, we measure the misalignment as a function of time. We 
define the inner disc as particles with r < 3 comoving kpc and we 
measure the halo properties in a shell between 0.2 < r/rvir < 0.5. 
We limit the shell to 0.5 r/r^ir to minimize the variations due to 
the newly infalling matter. In the left and centre panels of Figure [8] 
we show the orientation of the gas halo (solid lines) and the main 
disc (dashed lines) with respect to the static simulation frame as a 
function of time. In the right panel we show tp, the angle between 
the angular momentum vectors of the two components defined in 
equation[T] 

The two components enter the quiescent period of evolution 
misaligned by > 50°, following the last major merger at z ~ 2. 
The disc's orientation is essentially constant between 4-9 Gyr prior 
to the minor merger, which takes place at ~ 9.5 Gyr, but the halo 
orientation is slowly changing due to ongoing accretion. The minor 
merger provides an important clue as to why the inner disk and the 
outer gas halo are misaligned in the first place. The infalling satel- 
lite rapidly torques the disc, changing its orientation in the inertial 
frame and with respect to the hot gas halo. The gas halo, however, is 
unaffected by the merger because the infalling satellite is not mas- 
sive enough to have a considerable effect on the halo structure or 
angular momentum content. The misalignment between the differ- 
ent components is therefore perpetrated by such interactions that 
affect the inner and outer regions differently. 

Figure [3] shows that the warp ends at ~ 1 1 Gyr, yet Figure [8] 
shows that even z — the gas halo is misaligned with respect to 
the central star forming disc. If the warp is due to the misalignment, 
why does it disappear while the misalignment persists? Note that 
although the warp is long-lived, the gas never lingers there for more 
than 1-2 Gyr and it therefore must continually be supplied from the 
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Figure 8. Evolution of direction of angular momentum vectors of the main disc (r < 3 kpc) and the gas halo (0.2 < r/r^ir < 0.5). 9 and are angles 
between the angular momentum vectors and the z and x -axes of the simulation box, respectively. t/> is the angle between the angular momentum vectors of the 
two components. During the major merger at ^ ~ 2, the evolution of the spin axes of the two components is decoupled, resulting in an overall misalignment 
between the central disc and the rest of the gas halo. 
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Figure 9. Mass of infalling gas with T < 10^ K found in a 1 kpc-wide 
shell at 30 kpc from the centre. Assuming that this gas will reach the disc 
in a short time, we may use it as a proxy for an accretion rate onto the disc. 
The sharp drop at ~ 11 Gyr is coincident with the disappearance of the 
warp, therefore implying that the appearance of a large-scale warp requires 
ongoing accretion of cold gas. 

halo (see Figure|6l left panels). Thus, a gaseous warp may disappear 
when it ceases to be resupplied from the halo. 

To test this hypothesis, we evaluate how much gas is being 
supplied to the disc as a function of time. In Figure[9]we show the 
mass of cold gas falling toward the centre through a 1 kpc shell 
at 30 kpc (comoving) away from the disc. The shell is chosen to 
be far enough from the centre that it does not include any part of 
the disc, yet close enough that we can reasonably expect the cold 
gas in this shell to make it to the disc in a short time. We show the 
contribution of unshocked and shocked particles to this cold gas 
reservoir by dotted and dashed lines respectively; the remainder is 
made up of gas that enters the halo as a part of a merging event. 

Figure [9] shows that a significant drop in the amount of disc 
fuel coincides with the disappearance of the warp, confirming our 
expectations that the warp is kept alive by continual accretion. Un- 
shocked gas is most important for disc growth in this galaxy (see 
Figure 7 of Brooks et al. 2009), therefore it follows that the de- 
crease in the unshocked fraction has the largest effect on the warp. 
Further contributing to the disappearance of the warp is the growth 
of the central disc. As it grows, it begins to dominate regions where 
the warp had been prominent before. 



We note that satellite accretion events or perturbations do not 
necessarily result in large-scale warps. At ~ 12 Gyr, there is an 
influx of gas stripped from a satellite, shown by the excess cold gas 
which is not attributed to either shocked or unshocked components, 
but the perturbation fails to incite a warp. The cold gas contributed 
by the satellite by itself does not reach requisite densities to appear 
in the warp. 



5 PROPERTIES OF STARS IN THE MISALIGNED DISC 

Extended star formation in the misaligned gas component must 
leave a trace in the stellar populations at the present day. We can 
begin to investigate the importance of stars forming in the mis- 
aligned component by dividing stars based on whether they were 
formed interior or exterior to the warp radius at their time of for- 
mation. Determining the galactocentric radius of formation in the 
disc plane is not trivial in a cosmological simulation because nei- 
ther the disc's centre nor its orientation in the simulation frame are 
fixed. Further, the variations of the disc's orientation are tempo- 
rally poorly sampled because for practical reasons we only retain a 
limited number of outputs. Simply extrapolating the disc's position 
and orientation at each snapshot and the time of formation of stars 
proves insufficient. 

However, we do retain full information about the positions and 
times of star formation events for the entire simulation, regardless 
of when the outputs are produced. To proceed, we make the as- 
sumption that in the main galaxy, the vast majority of the star for- 
mation occurs in the disc and is symmetric about the disc's centre. 
With this assumption, and given a sufficient amount of star forma- 
tion, we now have knowledge of the disc's 3D orientation within the 
simulation volume, sampled at the star formation timescale (~ 1 
Myr in our recipe). We bin stars according to their time of forma- 
tion in 10 Myr bins and for each bin determine the disc plane. With 
this finely sampled orientation information, we then determine the 
precise radius of formation for each star in the plane of the disc at 
that time. 

In Figure [TO] we show the star formation history in the mis- 
aligned disc and main star-forming disc. We restrict our analysis to 
stars with disc-like kinematics at z=0, defined to be stars that satisfy 
Jz / Jcirc{E) > 0.8, where Jcirc[E) is the angular momentum of 
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Figure 10. Stai' formation rate as a function of time in the misaligned and 
main star-forming disc. In blue we show the star formation rate from cold 
gas accretion, in red from hot accretion, and in green from clumpy accre- 
tion. See text for definitions. 



a circular orbit at the particle's energy. Colours correspond to the 
accretion mode of gas, as in Figure[3] 

Figure [Tol shows that the unshocked gas is the dominant fuel 
for star formation during most of the evolution, in both the main 
and misaligned discs, which disappears almost completely at 11 
Gyr. Until the last few Gyr, the SFR in the misaligned disc is lower 
than in the main disc by factor ~ 2, although the total gas mass 
is higher (see Figure [S}- This is due to overall low gas densities 
in the outer misaligned component and correspondingly low star 
formation efficiency. 

The inefficient SFR in the misaligned disc implies that the en- 
richment of the gas in the misaligned disc will proceed much more 
slowly than in the main disc, resulting in a more metal-poor popu- 
lation. In addition, whatever enrichment takes place, the metallicity 
of the misaligned disc is continually diluted by the infall of pristine 
gas, further keeping the metallicity low. We show this difference 
explicitly in Figure [TT] where we plot the metallicity distribution 
function (MDF) for the stars forming in the main disc and those in 
the misaligned disc. The offset in the peaks of the distributions is 
~ 0.8 dex. The detailed differences between the two MDFs will 
vary depending on the properties of the misaligned disc. However, 
the short life-cycle of gas in the misaligned disc should be general 
and will therefore always result in low enrichment. 

The top panel of Figure [T2l shows the z=0 radial surface den- 
sity profiles of all of the star particles making up the disc (solid 
line), decomposed into those stars that formed in the misaligned 
disc (dashed) and those formed in the main star-forming disc (dot- 
ted). The absence of misaligned-disc stars in the centre is due to 
our definition of misaligned disc stars as having Rtirth > 3 kpc. 
Stars that formed in the misaligned disc completely dominate the 
present-day outer disc. This is natural if we consider that the mis- 
aligned disc component persists for 7 Gyr and is by definition al- 
ways beyond the main star-forming disc. 

Figure[T2]suggests that the stars forming in the warp dominate 
the present-day surface density at i? > 6 kpc, or roughly 2 scale- 
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Figure 11. Mass-weighted metallicity distribution function for the stars 
found in the main disc (solid line) and the misaligned disc (dashed line). 



lengths. This is partly due to the fact that during the lifespan of the 
warp, the main disc region extends only to 5 kpc (only ^ 2 scale 
lengths). The warp may therefore be dominant in a larger fraction 
of the disc in this model than in some local examples, where warps 
are seen to begin only at several disc scale lengths. Consequently, 
the amount of stars formed in the warp is likely too high in this 
model, and their contribution to the thick and extended stellar popu- 
lation may be unusually large. However, we do not have constraints 
from more distant systems, where it is possible that warps begin at 
smaller galactocentric radii. 

Stars formed in the warp dominate the stellar density at larger 
radii, and because their star formation history is rather extended, 
the mean age profile shows a slight upturn in the outer disc. This 
upturn in the age profile is similar to that in R08, which was caused 
instead by radial migration. In cosmological AMR simulations, 
[Sanchez-Blazg uez et al. (2009) also found that a change in the age 
profile was due to prolonged star formation in the outer disc. How- 
ever, as we discuss in §[8]below, we cannot use the present models 
to determine the relative importance of stars forming in the warp to 
stars that may have migrated there, because the discs presented here 
do not support the spiral structure required for migration. For com- 
parison, the migrated population at 15 kpc in R08 resulted in a sur- 
face density of ~ 10''Mq kpc^^, which is an order of magnitude 
higher than the amount of stars we find forming in the misaligned 
disc. We therefore stress that while stars may indeed form in the 
outer disc over extended periods of time, given a disc that produces 
realistic dynamical behavior, the migrated population should con- 
tribute substantially to the stellar density in the outermost parts of 
the disc. However, in lower mass spirals where we may expect spi- 
ral structure to be less important than in Milky Way-size discs, the 
migrated population may well be sub-dominant in the outer disc. 

The star formation in the warp is also problematic because 
much of it occurs near the gas density limit imposed by the sub-grid 
star-formation prescription. We must therefore allow for the possi- 
bility that star formation in the warp is somewhat overestimated in 
our model. However, large-scale warps such as the one described 
in this Paper are common in many cosmological simulations. Be- 
cause of their potential to dominate over the in-plane population at 
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Figure 12. Top: Surface density profile sliowing all disc stars (solid), stars 
formed in the misaligned component (dashed) and stars formed in the plane 
of the main disc (dotted). Stars formed in the warp dominate the outer parts 
of the disc. The absence of stars formed out of the plane at low radii is due 
to our method of selection. Bottom: Mean age profile of disc stars. 

a few scale-lengths, stars formed in the warp should be treated with 
caution, especially when analyzing the overall kinematic proper- 
ties of discs formed in cosmological simulations as they are likely 
to exacerbate the problematic high velocity dispersions (see §[8]l. 



6 THICK DISCS FROM OUTER DISC MISALIGNMENTS 

The dominance of stars that formed from the misaligned gas com- 
ponent in the outer parts of the disc in our simulation opens the 
interesting possibility that such stars could contribute substantially 
to the present-day out-of-plane population. Because the disc grows 
with time, what was the outer disc at z=l is a modes t distance from 
the centre at z=0. In Figure[T3]we show SUNRISE ( Ijonssonll20o3) 
edge-on I-band images of disc stars at z=0, split in the same way 
as in Figure [12] The left and right panels show stars formed in the 
main star-forming disc and in the misaligned component respec- 
tively. Although they have disc-like kinematics, the stars which 
formed out of the disc plane constitute a much thicker present-day 
population than the stars bom in the main disc. Though stars form- 
ing in the warp at any given time are in a tilted disc, their orbits 



precess resulting in a thickened population by present day rather 
than an inclined disc. We defer making detailed comparisons with 
observed thick discs to future work, but we point out that they could 
be significantly affected by long-lived warps. 

Given the tenuous nature of the misaligned component man- 
ifested in its low star formation efficiency, the amount of stars 
forming out of that gas is somewhat dependent on our choice of 
star formation parameters. The parameter set used here has been 
shown previously to yield excellent a greement with obser vations 
in terms of the Schmidt-Kennicutt law dStinson et al.ll2006h . mass- 
meta llicity relation terooks et alj|2007t) and the TuUy-Fisher rela- 
tion jGovernato et al.ll2009h . providing evidence that the SF recipe 
used with these parameters is remarkably successful in a globally- 
averaged sense. However, star formation in such low-density re- 
gions in the observed universe may occur in localized clum ps (e.g. 
iFereuson et al.ll998l : lThriker et al.l20ol:lDong et al.l2008h and it is 
not clear that at the present state of simulations we can satisfacto- 
rily resolve such processes for Milky Way-size galaxies. Recent at- 
tempts to resolve molecular clo ud masses in cosmological simula- 
tions yielded promising resu lts jAgertz et alj|2009| j Ceverino et al.l 
I2OO9L lOovemato et alJ[2oT^ . but due to computational cost only 
the simulations of dwarfs have been carried out to z = 0. To fully 
understand the impact of the misaligned component on the disc 
stellar populations, a full parameter study of the star formation pre- 
scription is required, using higher resolution simulations than those 
currently available to us, where star formation will be resolved in 
localized clumps. 



7 PRESENT-DAY OBSERVABLE CONSEQUENCES OF 
MISALIGNED ACCRETION 

In the previous sections we discussed the consequences of star for- 
mation in the warp for a system whose warp disappeared a few Gyr 
ago. However, many galaxies in the local universe show large-scale 
gas warps, so we now turn to a simulation whose warp is preserved 
to the present-day. The details of this model are described in §|2] We 
have verified that the physical reasons for the warp in this model are 
analogous to ones in our fiducial simulation described in the pre- 
ceding sections. Although we do not show it here, a Briggs figure 
for this model (like those shown in Fig.|2ll traces out a smooth lead- 
ing spiral. Below, we describe the observational predictions mani- 
festing in a unique stellar population signature from warps forming 
by misaligned gas infall. These predictions should make it possible 
to directly determine whether nearby large-scale warps are caused 
by such a mechanism and therefore confirm the presence of both, 
current cold accretion and the existence of a hot gas halo. 

In Figure [14] we show the projected stellar density of stars in 
this model at present day with overlaid HI column density contours. 
We split the stars into young (age < 200 Myr) and old (age > 3 
Gyr) and show the corresponding maps in left and middle panels re- 
spectively. Resolved star studies using the Hubble Space Telescope 
(HST) can dis cern between such stellar populations using CMD 
modeling (e.g. Ide Jong et al]|2007l) . The young stars are predomi- 
nantly found in the main disc and the HI warp, while the old stars 
dominate the more extended and largely spherical halo. 

In the right panel of Fig. [141 we show the age distributions of 
the two small boxes shown in the left and middle panels. The age 
distribution of the blue box, coincident with the warp, is dominated 
by recent star formation, though an older halo population is also 
present. The tilted and extended old population can be seen in the 
vicinity of the current warp in the middle panel, likely a relic of 
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Figure 13. Surface brightness maps of disc stars at z=G. The left panel shows stars formed in the main star forming disc, while the right panel shows stars 
formed in the misaligned component. Although the stars on the right were formed in-situ and have disc-like kinematics, due to their extraplanar origin they 
result in a thicker, lower surface brightness present-day population. 



star formation that had taken place in the misaligned component 
in the past. The off-warp box is completely dominated by an old 
population. 

Evidence for the processes described in this Paper may there- 
fore be found in similar large-scale warps, which have recently 
been associated with low levels of star form ation in the outer- 
most disc re gions of M83 and NGC5055 I Thilk er et all l2005l : 



ISancisietalJboOSh . Spectroscopic studies of ionized gas in the 
XUV disc of M83 indicate that the warped gas has low metallic- 
ity (0.1 - 0.2 Zq) consisterit^with enrichment from a stellar popu- 
lation younger than 1 Gyr 4Gil de Paz et al. 2007). A similar result 
has recently been re ported in a spectrosc opic study of a polar disc 
galaxy NGC4650A isoavone et alJ201(]t) . Studies of resolved stel- 
lar populations in and surrounding nearby warps could distinguish 
between different warp formation mechanisms and provide insight 
about the importance of on-going cold accretion and hot gas ha- 
los in massive spirals, both of which are inherently very difficult to 
observe directly. 



8 DISC HEATING AND DYNAMICS 

It is important to understand the relative importance of stars form- 
ing in the misaligned disc compared to other mechanisms that may 
also be responsible for substantial amounts of material in the outer 
disc. In particular, R08 showed that radial migration in a grow- 
ing galactic disc can efficiently populate the outer disc, making the 
outer disc a potentially valuable laboratory for studying such pro- 
cesses. Can we use the cosmological simulations to understand how 
secular evolution compares to the processes operating on a larger 
scale? Radial migration results from the interaction of star particles 
with transient spiral arms, f ormed through gravitational instabilities 
dSellwood & Binnevll200A . However, in the absence of a cooling 
mechanism, a single spiral will heat the di sc enough to stabilize 
it com pletely within a few rotation periods dSellwood & Carlberd 
Il984h . In the simulations of R08, the disc was cooled by the con- 
tinual accretion of gas and formation of kinematically cold stars, 
thus remaining at the cusp of instability and providing conditions 
appropriate for continual transient spiral triggering. R08 showed 
that even when we consider the change in radius for particles on 



mostly circular orbits we see that many particles are far from their 
birth radii by the end of the simulation. Thus the observed differ- 
ences between the galactocentric radii of stars at the end of the 
simulation and their birth locations are not only because of orbital 
heating, but rather due to significant changes in the guiding centers 
of the stars. 

We now try to determine whether the disc simulated in a fully 
cosmological setting described in the preceding sections is subject 
to the same radial migration mechanisms seen in the idealized sim- 
ulations of R08. Figure[T5lshows the conditional probability distri- 
bution of formation radii given the final radii for particles formed in 
the main star-forming disc. We only plot a kinematically-cool sub- 
set of the disc particles, defined by having Jz/ Jcirc{E) > 0.95, 
where Jcirc(E) is the maximum angular momentum in the disc 
at that particle's energy. This selection ensures that the spread 
around the 1 : 1 line in Figure [T5] is due to secular redistribution 
of star guiding centers rather than disc heating. Figure [15] shows 
that very little radial migration has taken place in the presented 
model, with the 75 per cent mass contour indicating at most 1 kpc 
difference between R final and Rjorm (in the simulations of R08, 
this difference is ~ 4 kpc at 8 kpc). This differs from the recently 
reported results studying radial migration in AMR simul ations 



reported results studying radial migration m AMK simul ations 
dSanchez-BlazQuez et al.ll2009l : iMartinez-Serrano et al.1 120090 . al- 



though we emphasize that care must be taken in cosmological sim- 
ulations to distinguish heating from migration. 

Given that the simulations of R08 were idealized models, 
it is important to examine the underlying reasons for the dif- 
ferent behavior found in the cosmological simulation. It is a 
known problem of cosmological simulations that the resulting discs 
are too kinematically hot, suppressing the formation of asymme- 
tries. In Figure 1161 we sh ow the Toomre Q stability parameter 
teinnev & Tremainell2008h 



Q 



3.36EG 



(5) 



for the disc of R08 and the MWl disc 5 Gyr ago and at present- 
day. We calculate Q using only the stellar component. The Q pa- 
rameter for the isolated dis c is between 2-3 (marginally stable). 
ISellwood & Carlberd (Il984l) found that in a model where the disc 
is continually cooled, the value of Q reached a quasi steady-state 
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Figure 15. Probability density of formation radius as a function of tinal 
radius for stars that formed in the main disc plane. Cell values are scaled 
such that an integral over the values in each column multiplied by the area of 
a cell equals unity. The white contours enclose 50 per cent and 75 per cent 
of the mass at each given R final- We find little evidence for significant 
radial redistribution of stars. 



of > 2, which allowed for continual formation of transient spi- 
rals, in qualitative agreement with the behavior seen in the R08 
isolated simulation. The disc of MWl, however, heats to signifi- 
cantly higher values at modest radii and most of the disc is at Q 
> 3 even at present day when the disc extends to > 10 kpc. Such a 
disc carmot be expected to support spiral structure required for ra- 
dial migration. Indeed, the amplitude of the m=2 mode for this disc 
is < 0.02 everywhere beyond the bar at z=0, a degree of az i muth al 
asymmetry appropriate for early-type discs IrIx & Zaritskvlll99^ . 
In the model of R08, the m=2 amplitude ranges from 0.02 to 0. 1 
at the end of the simulation. The total velocity dispersion {atot ~ 
+ cr^ + cFz) for a 1 kpc ring centered at 8 kpc exceeds 100 
km s^^ for the cosmological model, which is much higher than 
the measured value of 70 km s~^ fo r the oldest stars in the solar 
neighborhood l lHolmberg erai]|2009h . 

The current state of cosmological simulations available to us 
therefore does not allow us to address the relative importance of 




Figure 16. Toomre Q parameter as a function of radius for the MWl disc 
at two different epochs, compared with the Q measured in the isolated disc 
of R08. The cosmological disc is kinematically much hotter at most radii, 
which prevents the formation of spiral structure. 



radial migration compared to star formation in warps for the outer 
disc stellar populations. However, we feel justified in addressing 
the more general issues of angular momentum alignment and ac- 
cretion, since the relevant processes are seeded in the larger-scale 
cosmological dynamics. 



9 DISCUSSION AND CONCLUSIONS 

The high frequency o f warps in the nearby universe (e.g. 
iGarcia-Ruiz et al.l l2002h warrants efforts to understand the 
phenomenon in some detail. The range of warp morphologies 
is vast, though empirical rules do exist for the behavior of 
observed warps (Brisss 1990). In this Paper, we presented a 
possible formation mechanism for large-scale warping of disc 
gas in a Milky-Way size galaxy, similar in spirit to previous 
theories in that it relies on cosmological accretion of mate- 
rial. However, virtually all previous theories regarding warps 
postulate the warping of di sc mater ial after it had settled in 
the main disc plane (e.g. 



my, j m t... pim.>- yv-ft. Toomrd 1 19831: lO striker & Birme 

11989"; IPeb attista & Sellwood 1 1 9991: Ishen & Sellwood 



iinnevi 
20061: 



Weinberg & Blitzl l2006l : iDubinski & Chakrabartvl |2009|) . In this 
Paper, we have described a fundamentally different warping mech- 
anism, one which ties directly the morphological appearance of 
the gas disc to the angular momentum properties of the (invisible) 
gas halo. In this mechanism, the warp is a result of the hot gas 
halo torquing the infalling cold gas into a misalignment with the 
central disc. The central disk and the hot gas halo are misaligned 
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because they respond differently to accretion and substructure 
perturbations. Tfie warp is tfierefore at once an indication of 
ongoing cold gas accretion as well as the spin of the hot gas halo 
outside the central region. 

"Cold accretion" discussed here should not be confused with 
high- r edshift "cold streams " (e.g. IPekel et alJl2009al : IXgertz et al.l 
I2OO9I : ICeverino et alj|2009 '). which are a result of filaments pene- 
trating the virialized halo directly to the central disc region. The 
cold accretion at low redshift is instead in cold clouds and not fila- 
mentary in nature inside the virial radius. Nevertheless, the fact that 
some gas enters the halo at relatively low temperatures < ICK has 
large consequences for the disc-building process, as evidenced in 
the preceding sections. 

Not all galaxies show warps of quite the same magnitude 
as the one we have described in this Paper, and in particular 
the current amplitude of the Milky Way warp appears smaller. 
We restricted our discussion to the large warping specifically 
to focus on its cosmological implications. Mechanisms for gen- 
erating smaller vertical offsets in the disc plane appropriate for 
the current state of the Milk y Way warp h a ve been discussed 
by m any other authors (e.g. Toomre 1983; Ostriker & Binney" 
19891: IPebattista & Sellwood 1999; Shen & Sellw ood, ,2006 ; 
Weinberg & Blitzll2006l ; iDubinski & Chakrabartvll2009h . 

Although our own Galaxy does not show such a prominent 
warp, extragalactic systems of this sort abound in the nearby uni- 
verse. Much effort has in recent years been devoted to uncover- 
ing sources of disc fueling, such as isolated high-velocity clouds 
(HVCs) either around the Milky Way or other galaxies. The gen- 
eral conclusion is that observed sources of HI gas do not ac- 
count for s ufficient accretion rates to sustain current observed SFRs 
dSancisi e t al. 2008). We have shown that the large warp requires a 
constant supply of fresh cool gas (Figure |9j, and argued that the 
lifespan of individual parcels of gas in the outer disc is short. Thus, 
the presence of such a feature should be considered together with 
other indicators of accretion, such as discrete cold clouds, as evi- 
dence for ongoing cosmological gas infall. 

We have shown in § |7] that the warp formation process de- 
scribed in this Paper may be observable in nearby systems. Specif- 
ically, we have argued that the connection between XUV discs and 
HI warps is a direct result of the warping processes described in 
this Paper. Furthermore, the unique stellar population signature im- 
printed by this particular warping process should be discernible in 
inclined nearby systems where resolved star observations of outer 
discs are possible. 

We stress, that large-scale warps and their effects on stellar 
populations have largely gone unnoticed in recent work on cosmo- 
logical disc galaxy formation, although warps readily form in most 
such models. The amount of stars forming in the misaligned disc is 
likely model-dependent regardless of the simulation method used, 
but stars formed in this way can drastically alter the stellar popula- 
tion properties in the outer parts of the disc, including kinematics 
and metallicity (see §|5]and §[8j. Although we used a single simu- 
lation in this work to describe the physical process responsible for 
the warp, many cosmological simulations of isolated spiral galaxy 
formation yield a large-scale warp during some part of the galaxy's 
evolution. The details of the misalignment (longevity, degree of tilt) 
will differ due to different merger histories, but the general features 
outlined here should be valid for all systems of similar mass. 
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